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ABSTRACT 



Context. Star formation rate (SFR), metallicity and stellar mass are within the important parameters of star-forming (SF) galaxies that 

characterize their formation and evolution. They are known to be related to each other at low and high redshift in the mass-metallicity, 

mass-SFR, and metallicity-SFR relations. 

Aims. In this work we demonstrate the existence of a plane in the 3D space defined by the axes SFR [log(SFR)(MQyr"')], gas 

metallicity [12+log(0/H)], and stellar mass [log(Msta,/Mo)] of SF galaxies. 

Methods. We used SF galaxies from the "main galaxy sample" of the Sloan Digital Sky Survey-Data Release 7 (SDSS -DR7) in the 

redshift range 0.04 < z < 0. 1 and r-magnitudes between 14.5 and 17.77. Metallicities, SFRs, and stellar masses were taken from the 

Max-Planck-lnstitute for Astrophysics-John Hopkins University (MPA-JHU) emission line analysis database. 

Results. From a final sample of 32575 galaxies, we find for the first time a fundamental plane for field galaxies relating the SFR, gas 

metallicity, and stellar mass for SF galaxies in the local universe. One of the applications of this plane would be estimating stellar 

masses from SFR and metallicity. High redshift data from the literature at redshift ~0.85, 2.2, and 3.5, do not show evidence for 

evolution in this fundamental plane. 
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Relations between important properties of astrophysical objects 
often lead to the discovery of the so-called fundamental planes 
when three parameters are involved. The fundamental plane (FP) 
for e lliptical galaxies (Djorgovski & Davis 1987; Dressier et al. 
11987) . relates their luminosity, velocity dispersion (dynamics), 
and scale length (morphology). This FP represents an impor- 
tant tool to investigate the properties of early type and dwarf 
galaxies, to perform cosmological tests, and compute cosmolog- 
ical parameters. It is also an important diagnostic tool for galaxy 
evolution and mass-to-light (M/L) variations with redshift. 

Fundament al planes have a l so be en defined for glob- 
ular clusters (Mevlan & Heggie '1997) and galaxy clusters 
(ISchaeffer et al. 1993: Adami et al. 1998). The parameter space 
of globular clusters, elliptical galaxies, and galaxy clusters is 
properly described by a geometrical plane L oc R^cr^, where L is 
the optical luminosity of the system, R is a measure of the size 
scale, cr is the velocity dispersion of the system, and a and /3 are 
free parameters. The FP for globular clusters, elliptical galax- 
ies, and galaxy clusters have very similar slopes, which means 
that, accounting for differences in zero points, a single FP with 
a range of about nine orders of magnitude in luminosity c an be 
defined (Schaeffer et anil993Hiba^a & L6pez-Cru3l2009l) . 

The FP that we introduce here relates three fundamen- 
tal parameters: the SFR [log(SFR)(M0yr"')], gas metallic- 
ity [12H-log(0/H)], and stellar mass [log(Mstar/M0)] of field 
SF galaxies. All these variables have been related in the 
past b y the mass-metallicity (M - Z) relation ([Leq ueux et d] 
Il979l) . the mass-SFR relation dBrinchmann et al.i i2004l) . 



and the metal licity-SFR relation (iLara-Lopez et alj 120101 : 
iLo pez-SancheJ |201(^) . Also, some autho rs have studied the 
int er-dependence of those variables (e.g . iHoopes et al.1 120071: 
Ellison et al.ll2008t iMannucci et a l."2010'). In the present work, 
we propose the generalization of those relations defining a plane 
formed by a linear combination of two of those variables with 
respect to the third one. 

The M - Z relates the mass and metallicity of galaxies, with 
massive galaxies showing higher metallicities than less massive 
ones, and it has been well established for the local universe (z ~ 
0. 1) by the work of Tremonti et al. ( 2004) using SDSS data. The 
M - Z relation ha s also been studied at low redshifts z ~ 0.35 
JLara-Lopez et alj I2009a.b.) . at intermediate redshifts z ~ 0.7 



(e.g.. iRodrigues et aP 2008|). and at high red shift z ~ 2.2 and 
z ~ 3.5 (lErb et al.ll2006HMaiolino et al.t2008l respectively). 

The stellar mass of SF galaxies is also related to the 
SFR, i n the sense that more massive galaxies s how higher 
SFRs ("Brinchmann et al.' '2004^, ISalim et all l2005h . However, 
Brinchmann et al. (2004) emphasized that at log(Mstai/MQ) > 
10, the distribution of SFRs broadens significantly and the cor- 
relation b etween stellar mass and SFR breaks down. At higher 
redshifts, iNoeske et al.l (l2007h showed the existence of a "main 
sequence" (MS) for this relation over the redshift range 0.2 < 
z < 1.1, with the slope of the MS moving as a whole as z in- 
creases. 

The metallicity and SFR of SF galaxies are weakly corre- 
lated, as will be observed in Fig. 1. However, and despite of the 
high scatter, SFR increases with metalhcity (iLara-Lopez et al.l 
120 1 Ot lL6pez-Sanchezl20 1 Ol) . 
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This paper is structured as follows: Sect. 2 describes the data 
selection as well as the SFRs, metallicities and stellar masses 
estimations given by the Max-Planck-Institute for Astrophysics- 
John Hopkins University (MPA- JHU) groupQ and adopted in this 
work. In Sect. 3 we define the FP for field galaxies, and conclu- 
sions are given in Sect. 4. 



2. Data processing and sample selection 

Our study was c a rried out with galaxies from SDSS-DR7 
(lYork et al.1 120001: lAbazaiian etanl2009) . Data were taken 
with a 2.5 m tele scope located at Apache Point Observatory 
(iGunn et al.l l2006l) . further technical details can be found in 
Stoughton et al. (2002). 

We used the emission-line analysis of SDSS-DR7 galaxy 
spectra performed by the MPA-JHU group. From the full dataset, 
we only consider objects classified as galaxies in the "main 
galaxy sample" (Strauss et al.l l2002l) with apparent Petrosian r 
magnitude in the range 14.5 < r < 17.77 and redshift range 
0.04 < z < 0.1, that represent a complete sample in magni - 
tude and redshift (iKewlev & Emsonll2008t iKewlev et al.ll2006h . 
The lower limit in redshift ensures covering > 20% of the 
galaxy light, which is the minimum required to avoid domina- 
tion of th e spectrum by aperture e ffec ts (Kewlev et al. 2005). 
Follow ing IKewlev & Ellisoiil (l2008h and lKobulnickv & Zaritskvl 
(11999 *). for reliable metallicity estimates we selected galaxies 
with a signal-to-noise ratio (SNR) higher than 8 for the Ha, H/?, 
[Nil] .16584, [O ii] ^3727, [Oiii] ^5007, and [S ii] AA 6717, 6731 
lines. However, using less restrictive criteria does not affect the 
relations derived here, but only incre ase their dispersion. For a 
detailed analysis on the line SNR see lBrinchmann et al.l (l2004h . 
Fi nally, SF galaxies were se lected following the criterion given 
by Kauffmann et al. (2003a) for the BPT empirical diagnostic 
diagram log [Oiii] ^5007/H/3< 0.61/{log([Nii] /Hq')-0.05} + 1.3. 
From this final sample of 32575 galaxies, metallicities, stellar 
masses, and SFRs used in the present work were obtained by the 
MPA-JHU group following the methods described below. Since 
field galaxies are the dominant population of this sample, the FP 
presented here would be representative of field galaxies, and dif- 
ferent from the known FP of elliptical and clusters of galaxies. 

Metallicities were estimated statistically using Bayesian 
techniques according to Tremonti et al. (2004), based on simul- 
taneous fits of all the most prominent emission lines ([On], 
HjS, [Om], Ha, [Nii], [Sii]) using a mode l designed for the in- 
terpre tation of integrated galaxy spectra dCharlot & Longhettil 
1200 lb . Since the metallicities derived with this technique are 
discree tly sampled, th ey show small random off'sets (see for 
details pTremonti et al.l l2004). Any depen dence of SFR on the 
estimate d metallicity w ould be minor (JTremonti et al.l l2004t 
IBrinchm ann et al.l 1200 8). For this work, we selected galaxies 
with 12-i-log(0/H) > 8.4, corresponding to the upper branch 
of the R23. However, galaxies with 12H-log(0/H) < 8.4, corre- 
sponding to the lower b ranch of the R2 3 calibration, are poorly 
sampled, see Fig 6 of [ Tremonti et all (12004) and Fig. 1.1 of 
IKewlev & EUisonl (l2008l) . Therefore, to avoid a systematic dis- 
persion in the FP, and to work with an homogeneous sample, we 
selected galaxies with 12-i-log(0/H) > 8.4, which correspond to 
the ~99% of our SF sample. 

Total stellar masses were estimated from fits to the pho- 
tometry using the same modelling methodology as described in 
Kauffmann et al (2003), with only small differences with respect 
to previous data released. 



Finally, total SFRs for SF galaxies are derived directly from 
the emission lines, based on the careful modelling discussed in 
Brinchmann et al. ( 200 4), who modeled the emissi on lines in the 
galaxies following the Chariot & Longhettil (1200 ih prescription, 
obtaining a robust dust correction. Also, the metallicity depen- 
den ce of the Case B Ha/H;6 r atio is taken into account as well. 
The IBrinchmann et alJ (12004^) method offers a more robust SFR 
estimate than using, for example, a fixed conversio n factor be- 
tween Ha luminosity and SFR (e. g. lKennicuttlll998l) . 

The FP presented in this stu dy was initially identified 
by us using STA RLIGHT data dCid Femandes et all 120051: 
iMateus et alJ l2006) for the above described sample, but estimat- 
ing SFRs using the Ha luminosity and the Kennic ut t (. 1 998.) rela- 
tion, a nd metallicities following the calibration of lTremonti et al.l 
(2004*). Nevertheless, although the FP derived is the same, it is 
noteworthy that the plane has a lower scatter when using the ro- 
bust SFR and metallicity estimations derived by the MPA-JHU 
group. 



3. The fundamental plane 

As mentioned in the introduction, the SFR, stellar mass, and gas 
metallicity of SF galaxies are related to each other Their strong 
relation is evident when these data are plotted in a 3D space 
with ortogonal coordinate axes defined by these parameters. A 
careful inspection of this 3D representation, shows the existence 
of a plane (see Appendix A). 

The projection of galaxy data over any pair of the axes of 
this 3D space will reduce to the M - Z, metall icity-SFR, and 
stellar mass-SFR relations, as shown in Fig. 1 . iTremonti et al.l 
(2004) demonstrated that the metallicity of galaxies increases 
with the stellar mass in a relatively steep way from 10^ '^ to 10^°'^ 



Mq, bu t flattens above 10'°^ Mq (see Fig. la).lKewlev & Elhsonl 



http://www.mpa-garching.mpg.de/SDSS 



(l2008h recalibrated the M - Z relation of lTremonti et al.l (l2004l) 
using the completeness criteria of redshift and magnitud given 
in Sect. 2. However, as the Kewlev & Ellison (2008) fit (see 
Fig. 1.1 of that paper) departs for the low mass population, 
we fit a second order polynomial to the median mass in 
metallicity bins of the M - Z relation, 12-i-log(0/H)= aq + 
fli[log(M,t,,/Mo)]+fl2[log(M,tar/Mo)]2, withflo = -25.93923, 
fli = 6.39283, a2 = -0.29071, and (t = 0.26. Regarding the 
mass-SFR relation (see Fi g, lb), the SFR increase s with stel- 
lar mass up to ~ 10"^ M0 (IBrinchmann et alJl2004t) . while for 
higher mass values the scatter increases (see Fig. lb). We fitted 
a line to the median mass in log(SFR) bins of the mass-SFR re- 
lation, log(SFR)= flo-Hfli[log(Mstar/M0)], with ao = -12.50704, 
fli = 1.27909, and cr - 0.27. Note that the a given is that 
of the horizontal axes [log(Mstai7M0)] in both cases. However, 
for 12H-log(0/H) < 8.9, the SFR is not strongly correlated with 
metallicity (see Fig Ic), whereas for higher metallicity values, 
the SFR increases rapidly. Therefore, since the metallicity-SFR 
relation does not correlate at all, this leads us towards a possible 
solution: a linear combination of the metallicity and SFR as a 
function of the stellar mass. 

We generated least-squares fits of one and two order polyno- 
mials combining two of the variables as function of the third one, 
trying all the possible combinations. As expected, the best one 
found is a linear combination of the SFR and metallicity yield- 
ing the stellar mass (see Fig. 2). The FP for field SF galaxies is 
then defined by: 

log(M,„,/Mo) = a [12 + log(0/H)] + /? [log(SFR) (Moyr-')] + 7 (1) 
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Fig. 1. In panels a, b, and c, the M - Z, stellar mass-SFR, and metallicity-SFR relations are respectively shown. The density plots 
are represented in bins of 0.02 dex. Red circle s represent the median mass in bins of 0.05 dex in 12+log(0/H) and of 0.13 dex in 
log(SFR). Blue empty diamonds represent the' Erb etal.1 ( 12006 ) data at z~2.2, the bars of their data show the metallicity error and 
the stellar mas s range of each bin, while the SFR errors are of the order of ~ 1.3 dex. Red filled diamons represent the sample of 
iMaiolino et al.l (l2008l) at z~3.5, the bars of th eir data show the metallicity and mass errors. The black filled square represents the 
median of the data of 'Rodrigues et al.l (l2008h at z~0.85 with its respective error bars. Solid lines in panels (a) and (b) represent a 
two and one order polynomial fits to the median values, respectively. 
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Fig. 2. FP for field SF galaxies. The stellar mass, in units of solar 
masses, is presented as a linear combination of 12+log(0/H) and 
log(SFR). Symbols fo llow the same code as in Fig 1 . Since the 
IMaiolino et al.l (l2008l) data do not have SFRs errors, the vertical 
bars of their data show only the error in metallicity of the linear 
combination [a A(12+log(0/H))]. 



with a = 1.122 (±0.008), /3 = 0.474 (+0.004), y = 
-0.097 (+0.077), and cr = 0.16. The sigma given is that of the 
vertical axis of Fig. 2. 

One of the uses of this FP would be to estimate stellar masses 
from the metallicity and SFR of emission line galaxies. In anal- 
ogy to the FP of elliptical galaxies, it is expected that this FP 
reduces the scatter of the stellar mass found throught the M - Z 
or mass-SFR relations. To compare the mass (Mstai/Mo) width 
of this FP with that of the M - Z and mass-SFR relations, we 
estimated the 90th and ~68th mass percentile around the fit to 
each relation, obtaining a mass width for the 90th percentile of 
~ 0.57, ~ 0.87, and ~ 0.92 dex for the FR M - Z, and mass- 
SFR, respectively. Using a ~68th mass percentile, we found a 
mass range of ~ 0.32, ~ 0.52, and~ 0.53 dex for the same rela- 
tions. As observed, the scatter in the mass estimates through the 
FP is ~ 0.3 and 0.2 dex lower than that of the M - Z relation for 
the 90th and ~68th mass percentiles, respectively, which means 
that the mass dispersion using the FP is reduced ~38% for a 68th 
percentil. 



The largest expected errors of the stellar mass estimates that 
would be obtained using this FP, without considering SFR and 
metallicity errors, would be of +0.28 dex for a 90th mass per- 
centile, and of +0.16 dex for a ~68th mass percentile (Icr error). 
Therefore, as the FP is a well behaved linear relation, and its dis- 
persion in mass is lower than that of the M - Z and mass-SFR 
relations, it proves to be a useful tool for deriving stellar masses. 

In order to identify possible signs of evol ution of the 
FP, we used the data of iRodrig ues et al. (2008|) at z -0.85, 
lErb et all (l2006h at z -2.2, and of MaioHno e t'al] (12008.) at z 
~3.5, as shown in Fig. 1 and 2. The Erb et al] (l2006h metal- 
licity estimates, using the N2 method and the calibration of 
|Pettini & Page l (2004), and the Maiolin o et al. (2008), using the 
calibrations o fiKewlev & Dopital ( l2002|). were recalibrated to th e 
R23 method using the conversions of iKewlev & Ellisod (l 2008h . 
The lErbet all (2006) data use the Bruzual & Chariot (2003D pop- 
ulation synthesis models and a Chabrier (2003) IMF. According 
to Bru zual & Chario t (2003), t he Chabrier (2003) and KroupJ 
(120011) IMF, used bv lErbetall (12006) and in this study, respec- 
tively, yield practic ally identical M/L ra tios. However, since 
the ste llar masses of IMaiolino et al.l (12008) andlRodri^ues et al.l 
(12008 ^) are estimated using a Salpeter (1955) IMF, we cor- 
rect them, as indicated by Maiolino et al. (2008) , dividing thek 
masses by a factor of 1 . 17 to make these data consistent with our 
assumed IMF. 

In Fig. la the w ell known evolu t ion of the M-Z relat ion can 
be appreciated (e.g. lErb et aill2006l: JMaiolino et all2008l) . while 
in Fig. lb an important evolution of the log(SFR) of ~ 1.5 dex 
is observed for an intermediate mass of lO'^'M©, as well as a 
flattening in the sha pe of the mass-SFR relation, already noticed 
by Lara- Lopez et al. (2010). Th e evolution of the M-Z relation, 
as argued bv lLiang et al.l (l2006h is due to a decrease of the metal 
content in galaxies rather than an increase of their stellar mass. 
Also, the evolution observed in the mass-SFR relation is due to 
an increase of SFR in galaxies. Therefore, the metallicity-SFR 
relation shows evolution in both axes, as shown in Fig. Ic. 

We apply Eq. 1 to these high redshift data. As observed in 
Fig. 2, the Erb et al. (2006) data could show an evolution in the 
slope. However, the uncertainties of the measured parameters are 
too large. Moreover, if this change in slope were real, it would 
be also evident in the data of Maiolino et al. (2008), and it is 
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not. Therefore, we conclude that current high redshift data do 
not suggest an evident evolution of the FP. 

It is noteworthy that the main evolution in the M - Z and 
mass-SFR relations is driven by the metallicity and SFR, respec- 
tively, rather than the stellar mass. Therefore, we would not ex- 
pect a mass evolution in the mass proyection of the FP, as can be 
appreciated in Fig. 2. This lack of evolution in our relation could 
be explained by the metallicity and SFR evolving in opposite di- 
recctions, which means that high redshift galaxies would have 
lower metallicity values but higher SFRs compared to the local 
sample. For example, those diferences could be, for an interme- 
diate mass of -lO'O'Mo, of 2.12 dex in log(SFR), and of -0.7 
dex in 12+log(0/H). This means that the high SFRs at higher 
redshifts would be compensated with their lower metallicities 
when the a and/? coefficients of Eq. 1. are taken into account. 

Stellar mass is the physical fundamental parameter driving 
the SFR and metallicity of star forming galaxies. Both, SFR 
and metallicity increase with mass, as shown in the M - Z and 
mass-SFR relations of Fig. 1. A possible e xplanation is given 
by the w ell known effect of downsizing (e.g. 'Cowie et alj|1996t 
iGavazzi & Scodeggio 1996), in which the less massive galaxies 
form their stars later and on longer time scales than more mas- 
sive systems. This implies lower metallicities and higher spe- 
cific SFRs for low mass galaxies. Since star formation rate his- 
tory drives the metal enrichment, downsizing correlates both pa- 
rameters with mass. Then, their linear combination would relate 
current star formation rate with its past history, increasing the 
accuracy of the mass determination by reducing the scatter of 
the relation. 

Usually, the stellar mass of galaxies is estimated through the 
z-band magnitude, the spectral indices Dn(40 00) and H^^, and 
assum ing an IMF, such as the method used by iKauffmann et alj 
(l200 3b). It is also possible to use sophisticated codes such as 
STARLIGHT (Cid Fernandes et al. 2005; Mateus et al. 2006), 
which fits an observed spectrum with a combination of sim- 
ple stellar populations (SSPs) fro m the evolutionary synthe- 
sis mod els of Bruzua l& CharloJ yOOS), computed using a 
IChabrier (2003) IMF, and "Padova 1994" evolutionary tracks 
(Girardi et al. 1996). Moreover, the obtained masses must be 
corrected for aperture effects based on the differences between 
the total galaxy magnitude in the r band, and the magnitude in- 
side the fiber, assuming that the mass-to-light ratio does not de- 
pend on the radius (see Mateus et al. 2006 for details). Therefore, 
the use of the FP presented here allows estimating stellar masses 
in an easier and reliable way compared with the existent meth- 
ods. 

Fin ally, in an independent and parallel studv. lMannucci et alj 
(1201 Ol) fit a 2D surface instead of a plane to the same variables, 
which also reduces the metallicity dispersion. However, one of 
the most important differences between both studies is that we 
are using a complete sample in redshift (z<0.1) to avoid sys- 
tematics effects (e.g. Lara-Lopez et al. 2010). Additionally, the 
plane fitted by us also account for the data at higher redshifts, 
sustaining the apparent evidence of no-evolution laid down in 
this section, as explained above, which makes it useful for esti- 
mating stellar masses in field galaxies even up to z~3.5. 



4. Conclusions 

We have demonstrated the existence of a FP for field SF galax- 
ies in the 3D space formed by the orthogonal coordinate axes 
log(Msta,-/ Mo), log(SFR)(M0yr-i), and 12+log(0/H); three of 
the fundamental parameters of galaxies. All those variables have 



been related previously in pairs as with the M - Z, metallicity- 
SFR, and mass-SFR relations, but this is the first time that the 
correlation for all of them has been quantified. 

The FP presented here allows estimating the stellar mass 
[log(Mstai-/ Mq)] of field galaxies as a linear combination of 
12-i-log(0/H) and log(SFR)(M0yr"'). The scatter in the mass es- 
timates using the FP ( 1 cr error of 0. 1 6) is lower that that obtained 
through the M - Z and mass-SFR relations. 

The FP introduced here would be useful for deriving masses 
in spectroscopic surveys where the SFR and metallicity are 
estimated for emission line galaxies, for e xample, using th e 
Wa luminosity to estimate the SFR (e.g. iKennicutl Il998h . 
and any of the metallicity method s in the literature, su ch as 
the R;^ (lPageletalJll979h or N2 dDenicolo et all |2002|) . see 
iKewlev & Ellison! (l2008h for a review. However, since this study 
has been carried out using emission line galaxies, this FP will be 
useful only when both, SFR and metallicity of galaxies can be 
estimated. 

Within the errors, there is no evidence of an evolution of the 
local FP when applied to high redshift samples. Which means 
that it could be useful even at high redshifts, where measuring 
the continuum and absorption lines for fitting models would be 
more difficult and time consuming. 

Then, we propose the use of this FP as an alternative tool to 
the existing methods to determine the stellar mass of galaxies at 
low and high redshifts. 
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Fig. A.l. Tridimensional representation of the FP according to the fit given in Eq. 1, the color scale of the plane is related with the 
stellar mass. Light colour symbols are located above the plane, while the dark ones are below it. Dots represent the SDSS local 
sample described in Sect. 2, diamonds correspond to the Erb et al (2006) data at z~2.2, and circles symbolise the data at z~3.5 from 
Maiolino et al. (2008). The yellow triangle represents the data of Rodrigues et al. (2008) at z~0.85. The iso-density contours are 
subject to a projection of the SDSS data cloud on the metallicity-SFR plane, and they are given only as a visual aid. 



